Séparabilité de l'équation de Laplace en 3 et 4 dimensions 
pour un système de coordonnées non orthogonales entre elles 


Séparabilité d'un système de coordonnées sphéro-coniques non-orthogonales à 3 dimensions 


Le système de coordonnées sphéro-coniques possède la paramétrisation suivante : 


z2>0 y=vp z<0 u’ +v?’ <1 v= 


Ce n'est pas un système de coordonnées orthogonales. Une simple illustration de deux iso-surfaces 
en croisement permet de s'en rendre compte : 


Les paramètres u,v sont donc censés représentés une position quelconque sur la sphère de rayon p, 
moyennant la duplication des coordonnées sur les deux hémisphères. 


Dans ce type de paramétrisation, le plan équatorial z=0 est représenté par la relation u°+v°=1 et la 
variation du rayon p. Le plan équatorial est donc représenté par une suite de cercle concentrique 
p=Cste sur lequel se vérifie la contrainte u?+v/=1 .Sur ce même plan équatorial l'origine sera 
indifféremment représentée par toute valeur de u,v tant que p=0. 


La valeur u=1, contraint v=0 et correspond à l'axe des x. La valeur v=1, contraint u=0 et correspond 
à l'axe des y. L'iso-surface u=Cste correspond à un cône centré sur l'axe des x qui s'évase lorsque u 
décroît vers 0 pour devenir le plan (yz) . L'iso-surface v=Cste correspond à un cône centré sur l'axe 
des y qui s'évase lorsque v décroît vers 0 pour devenir le plan (xz) . 


Les éléments métriques de ce système de coordonnées vont nous permettre de représenter la 
forme du Laplacien, qui s'avère être le « Laplacien de Beltrami », formalisme utile pour se placer 
dans le cadre de coordonnées curvilignes générales (voir plus loin dans le document). La matrice 
Jacobienne s'écrit sous la forme : 


P 0 u 
J= 0 P v 


P P Me. 
u v l-u -v 
\1-u? -v \1-u°? -v 


Le tenseur métrique contravariant est le suivant : 


(i-ve /-u-v) (up? }/(1-u° -v°) 0 
g&=J J= (uv }/(1-u° -v°) (i-u)? -u -v?) 0 
0 0 1 


Le tenseur métrique covariant est le suivant (matrice inverse du tenseur métrique contravariant) : 


Le déterminant du tenseur métrique contravariant est particulièrement simple : 


4 2 
det(g,)= < zi Idet(g, | F P 


l-u° -v V1-u° —v 


L'équation de Laplace-Beltrami s'écrit alors : 


T= a ak Jdet(G) g'o re 


ar- Re l-u? ôT uv a "| uv OT. 1-v° z) 1 ð (2) 


p° Ou Vl-u°-vy? u J1-u?-v? ôv ðv 1-1 -v ôu siey ôv A y op op 
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Ou 1-1? -w Ou N1-u? -v ôv E 1-1? -vy Ou “iayy ðv êp Ôp 


J'exposerais plus loin une extension des polynômes de Legendre à deux dimensions permettant la 
séparation de l'équation de Laplace précédente. Mais avant cela introduisons les dès à présent. Les 
fonctions R°mn(x,y) sont définies comme suit : 


1 
-x-y E EE suis 
(a x+b y—1) + la? vo g y) 2 2.4 b" Ran (0y) 


1 
U? (x, y)=(l-x -y JR (57) 


Dont les deux équations aux dérivées partielles et la sommation de ces dernières sont les 
suivantes : 
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En revenant à l'équation de Laplace précédente, en développant les termes en u et v, il apparaît les 
termes de l'équation aux dérivées partielles de la fonction R°mn(X,y) : 


per ô l-u? ƏT uv ôT) ð uv oT 1-v° OT 
Ou 1-u?-v? Ou  J1-u?-v? Ôv ðv Vl-u° -y ôu E ôv 
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Ainsi si l'on sépare la solution recherchée de l'équation de Laplace sous la forme T= R, (u,v) Olp) 
alors la fonction Q est solution de l'équation différentielle ordinaire : 


, 


d (tel) Druide) 


dp dp 


Remarque : le terme (m+n+1)(m+n +2) équivaut à l'introduction de la variable de séparation dans 
l'équation de Laplace. 


Cette équation radiale possède la solution générale suivante : 


Q(p)= p° = a(a +1)=(m + n+1)(m+ n + D fe SRNE LL ofo)= ap" + Bot) 
a={(m+n+2) 


On peut donc proposer comme solution séparable dans ce système de coordonnées sphéro-conique 
la forme suivante : 


T{u,v,p)= E S 4,10" +B 


m=0 n=0 


aa ng (u, v) 


Ce type de développement est-il exploitable dans le cadre d'un problème aux limites, par exemple 
le plus « simple » comme celui de Dirichlet. La réponse réside dans le fait qu'au moins toute 
fonction f(u,v) ayant des propriétés suffisantes de régularité, ainsi que ses dérivées premières et 
étant bornées, ainsi que ces dérivées premières, dans le domaine u2+v2<=1 peut se représenter 
sous la forme d'un tel développement. Cela veut simplement dire que le système de fonction R est 
complet dans cet ensemble fonctionnel. 


Ce « Théorème » de représentation a été exposé dans deux ouvrages sur les fonctions hyper- 
sphérique et hypergéométriques à plusieurs variables, la thèse de 1914 de J.Kampé-de-Fériet « Sur 
les fonctions hypersphériques » chapitre IV « Développement d'une fonction arbitraire en série de 
polynômes zonaux », page 69 et celui de 1926 de P.Appell et J.Kampé-de-Fériet, « Fonctions 
hypergéométriques et hypersphériques, polynomes d'Hermite » chapitre V « Développement d'une 
fonction arbitraire en série de polynômes et fonctions V`mı.„mn(Xu..,Xn), U`mi..mn(Xu.-Xn) et 
R°m1...mNX 2e XN). 


Voici quelques détails introductifs de ces polynômes et du théorème de représentation. Les 
fonctions V mi.mn(Xu--,Xn), Umi.. m(Xu-- Xn) et R°m1.mnlXz.,Xn). Sont définis par leur fonctions 
génératrices suivantes : 


S mı =+0 My =+0 
2 2 2 2 2115 m; my 
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s 2e 2 2P rys ( ) (x, eT Xy 2 Es 3 m My DS ( ) 
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Soit sous forme de dérivations de la formule de Taylor appliquées aux fonctions génératrices : 
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Les fonctions et polynômes Vi..mvlXz..,XN), Umi. m(Xu. Xn) et R°m1..mnlX1.…,Xn) Sont également 
définis par des dérivées partielles suivantes : 


2 2 = 2 2 
Xy =l- x Xy Ey =l+é +..+6, U=M ++ Mr 


a Sn X% Xy 
m nE Te Èz P eE 
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De par leur construction V..mn{X1..,Xn), Um1.mvlX1..,Xn) sont des polynômes et R°,:1..mn{X1..,XN) 
une fonction qui s'annule sur l'hypersphère unité pour s>1 : 


s-1 


RÈ a (xpa xy)= Xy 2U 


mi My 


(xixu )=0 4 X, =06 x Este =l 


s 
m.. My 


Le théorème de représentation stipule que toute fonction F de plusieurs variables dans le domaine 
intérieur de l'hypersphère à N dimensions, ayant des conditions suffisamment « convenables » peut 
se représenter sous la forme d'un développement en polynômes comme suit : 


M=+0  My=+0 


F(x, Xy) = >. = Aam m (xis Xy) (x Xy )/ Xy =l- Xy >0 
0 


m=0 my= 


Les conditions suffisamment « convenables » portent sur l'existence d'un borne supérieure pour 
OHN F(a p.a Xy) 
da" dx, 
convergence du développement. Il se peut que le théorème de représentation puisse s'appliquer à 

un espace fonctionnel plus large, mais ce n'est pas étudié dans ces ouvrages. 


toutes les dérivées partielles de la fonction F : 


<H et permettent d'assurer la 


Les coefficients du développement se calculent à l'aide des intégrales suivantes : 


L=M +.+My Symbole de Pochhammer (s, u ) = 


T(s+u) (1, jar 


F(u) 


y (=) 
2m? 2 (s, u) K sa s 
À = |(X F3: 
on rf N tect) (ym (l,m) "7" K n)? Pro n) Unm 
2, 


Gin) dx,.…dxy 


Cette formule est justifiée par le fait que les deux jeux de polynômes V`mı.„mn(Xu..,Xn) et 


U°1..mvlX2..,XNn) s'avèrent deux à deux orthogonaux suivant le produit scalaire induit par la formule 
suivante : 


N s-l 
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Ce théorème de représentation s'applique par évidence au développement par les polynômes 
U°1..mvlX2..,Xn) car il suffit d'écrire : 


My=H0  My=+0 


F(x Xy) = È z È Baom Unom Gelée) ER )/ Xy =1-x n- Xy >0 


m,=0 my=0 
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Et également aux développements avec les fonctions R°mı..mn(Xu--,Xn) : 
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L'application à deux dimensions pour un développement de la forme : T(u,v)= 5 5 A R lur) 
m=0 n=0 ý 
Tlu, v) = € FaR, (u,v) 
donne donC: -2 e N=2 > FE 
À  mni{2m + 2n +3] f dudv T(u,v) y? (u,v) 
| 2a(m+n+1} Pa | 


Application à la construction du problème intérieur de Dirichlet de l'équation de Laplace sur la 
sphère 


x=up x=up 
Selon le système de coordonnées : |! =y p z220 y=vp z<0 
l-u’ -v° p z=-Vl-u -v p 


étant donné que la paramétrisation est double selon les z positifs et négatifs, il convient peut-être 


de scinder le problème aux limites de Dirichlet : AREO E ox = fe) en deux sous- 
problèmes selon le principe de superposition. Pour chacun de ces sous-problèmes, lorsque les 
conditions aux limites sont symétriques en z, cela signifie que la dérivée première de la fonction 
selon z doit s'annuler sur le plan équatorial 


Or la dérivée première selon z s'écrit : 


ÔT(x, y,z)_ ôu OT(x,y,z) Ov OT(x,y,z) Op ƏT(x,y,z) _ Iu? u OT(x, y,z) v ôT(x,y,z) OT(x. y,2) 
&z Oz Ôôu & a & Op p ĉu e a æ 


m=+0 n=+00 m+n+l 
MT) S 4.(£ ) R2 (uv) 


Avec un développement de la forme m=0 n=0 , le troisième terme de la 
dérivée i-u? -y? aT yz) o. Par contre le terme Jf, Flu OT(x, y,z) l gau lui ne 
op Ou ôv 


s'annule pas forcément du fait de la dérivée première du terme hga On écrit plus en détail 


terme à terme : 


Jida Baler, , Ball, EE Ph a! l-u’ -v U? „lu, DMG l-u?’ -v U? „(u,v À 


u ðv Ou Ov 


Pour 1-u° -v =0= V1-u° FL EE (u, T En kr (u, Dé tune -Un (u, v) 


ôu m,n ôv m,n 


Donc la condition pour que le gradient en z de la solution s'annule est que le développement en 
polynômes U°mn(u,v) s'annule de même : 


m=+%0 n=+00 m+n+l 
> > and2) U; (u, v)= 0 Vp et u’ +v’ =] 


m,n 
m=0 n=0 R 


Comme cette condition doit être vraie quelque soit le rayon p, on en déduit que la sous-sommation 
sur un même degré des polynômes en U’mn(u,v) doit être nulle, mais aussi qu'elle doit se vérifier 
quelque soit le jeu de coefficient, alors: Ym,n U? „(u,v)=0 pour u?+v?=1: Mais l'on peut 


démontrer que la forme des polynômes U°mn(u,v) est particulièrement simple : à une constante 
près il s'agit du produit d'un monôme en u et d'un monôme en v : 


Vm,n U? „lu, v)ocu”w" pour u +y =1 


qui ne peut donc être nulle sauf aux valeurs particulières u=0,1 soit v=1,0. En d'autre termes le 


m=+% n=+00 m+n+l 
> >. MA RÈ „(u, v) 


développement  "-0 »-0 R n'est pas adapté pour représenter une solution 
harmonique aux conditions limites symétriques en z. 


Lorsque le problème est anti-symétrique alors la fonction doit s'annuler sur ce même plan 


m=+0 n=+ m+n+l 
Tuvo) È TA (E) Rale) 


équatorial, or avec un développement de la forme: m=0 n=0 , la 


Vm,n RÈ „(u, v)=0 pour u’+vy’=1 


solution s'annule de fait car . Alors la solution du problème de 


Dirichlet : 
T(u,v, P) r = f (u,v) 


AT =0 
Fee =T(u,v,p)} =0 


u?+v?=1 


Est de la forme : 


HERO NES m+n+l 
Tu, v,p)= DNA, (£ ) R? (u,v) 


m=0 n=0 
m!n!(2m+2n+3) 
27z(m+n+1} 


m | dudv fluv) Vi, (uv) 
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Si maintenant la fonction limite f{u,v) sur l'hémisphère supérieure est pris comme une constante 
égale à 1 alors certaines intégrales s'annulent de facto du fait des propriétés de parité des 
polynômes Vhn{u,v). Ces derniers suive en u et v la parité de l'indice correspondant : 


Vg (eu, ev) = ("72 (uv) | dudv Vè,(u,v)=0 pour m ou n impair 


m,n m,n 
u?’ +u? <l 


La sommation ne se développe donc que sur les seuls indices m,n pairs, soit : 


M=+0 n=+%0 2m+2n+1 
np D È AR) Raate) 


m=0 n=0 R 


_ (2m) (2n)(4m+4n+3) 


u,v)=1= 4 
A) M 2rQm+2n+I} 


Í dudv F2 (u, v) 


u? +u’ <1 


À finir bien comprendre ce point, je suis rouillé ! Pourtant ce doit être simple !! et surtout tester 
numériquement 


Séparabilité d'un système de coordonnées axi-symétrique orthogonale à 3 dimensions 


Pour un système de coordonnées supposées orthogonales de type axi-symétrique autour de z (c'est 
à dire que toutes les iso-surfaces du système sont symétriques autour de l'axe z), la 
paramétrisation du système à 3 dimensions peut être représentée sous la forme : 


p (2) (me 


à ôp ôp 


x= F(p,0)Cos(p) 2 2 
»= F(p,0 )Sin(p) > ds? = hdp’ +h?do? +h?dp° 4h? = (122) leoh 
z= D(p,0) 

h, =(F(p,0)) 


Lorsque le système est orthogonal, Le laplacien et les éléments métrique ont la forme suivante : 


g = hhk? VE = ho SAT 1 (Er) l e(a), 1 (Ear) o 


prere | Je ol k õp) Jgôo\ k ôo) Jgõp\ hè õp 
h = 0F(p,0) Z 0D(p,0) | 
4 ôp ôp 


| See EE) _ à [rte r), ð (Fipo) r), hhe OT o 


ðc h, õp} ôo h, ôo) Flo,o)ôg? 
,0 


h, = 
h, =F(p 


Avec une solution séparable de la forme l'équation de Laplace restreintes aux coordonnées p,o 
devient : 


ô h, ôðU\ à h, OU) (m+1) 
T =U C l —| F <= — |+ —| F P h h U =0 
(o.oo) = Ulpio coslas o) = È] Plo. ESE), 2 Fipo) LE ut 


La condition d'orthogonalité est alors vérifiée lorsque les deux fonctions F et ® vérifient la relation 
suivante : 


oF(p,o)0F(p,0) 0D(p,o) 0D(p,0) 


=0 
Ôp 06 Ôp 06 


Extension à la séparabilité d'un système de coordonnées à 4 dimensions non-orthogonales 


Avant de rentrer dans le vif du sujet, faisons quelques rappels utiles de géométrie différentielle 
uniquement appliqués aux changement de variable conduisant aux systèmes de coordonnées 
curvilignes. Un simple exemple à 3 dimensions suffit pour comprendre. Nous allons d'abord rester 
dans le cadre d'un changement de variable de forme la plus générale possible et introduire les 
éléments différentiels utiles pour le manipuler. 


x= filea) q' = g'(x,,x3,x3) 

Xz = hleo) q sie) 

x, = f,(q',4°,q°) q = g° (xxx) 
On va supposer que toutes ces fonctions sont continues et dérivables dans le domaine des variables 
des deux jeux de variables. Les iso-surfaces du système de coordonnées curvilignes 
q'=Cste q°=Cste 4 =Cste s'appellent des surfaces de coordonnées, les courbes de 
coordonnées sont formées par les intersections de ces surfaces de coordonnées. Lorsque les 
coordonnées curvilignes sont dîtes orthogonales, alors les vecteurs tangents aux courbes de 
coordonnées, qui sont en intersection, sont orthogonaux entre eux. 


On peut construire une base de vecteurs indépendants dans le système de coordonnées qi par les 
dérivée première d'un point x dans l'espace : 


ôx, ôx, 0X 
ni Eee 
ôq ôq ĉq 
i ôx, Ox, x 
re) , [5 aa) 
ôx, ôx, àx, 
ce 
ôq” ôq” ôq 


Cette base est appelée base contravariante, elle peut être rendue normée par division par la 
norme. Cette base n'est pas à priori orthogonale. La base des vecteurs construite comme suit : 


p (ee a o 

CR êx, 

; ôq’ ôq’ ôq? 
Fa x 2. 3 k= Zie a, 
È k - c°) (2 ox, -0X3 
p [d ot or 

ax âx, Ôx, 


est appelée base covariante. On peut démontrer que cette dernière est orthogonale deux à deux 
avec la base contravariante pour les indices inégaux : 


Vi, j h.h, = ô; ô; symbole de Kronecker =1 sii=j =0sii+ j 


Un déplacement infinitésimal dans l'espace se mesure par la distance dx : 


dx = dq’ h, > ds? =|dx} = dg'dg'hh, 


Ici j'introduis la notation d'Einstein, pour laquelle lorsque les indices i sont communs à l'expression 
en position covariante (indice supérieur) et contravariante (indice inférieur) alors la sommation sur 
toutes les valeurs est sous-entendue. Il s'agit lorsque les formes différentielles n'implique qu'un ou 
deux jeux d'indice tout simplement d'un produit scalaire ou d'un produit matriciel. 


La matrice formée en colonne par la base contravariante de vecteurs est appelée la matrice 
Jacobienne. L'inverse de cette matrice est constituée en colonne par la base covariante de 
vecteurs : 


a mi OE A 
ôq! ôq? êg ôx, âx, ôx; 
2 2 2 
ee el pate 22 JJ’ =1=h h, =ô' 
ôq Ôq ôq ôx, âx, 0x, | 
Chr 0 og © ôq 
ôd af êr 6x, ôx, ôx, 


La matrice jacobienne exprime la transformation linéaire locale des vecteurs de base du système de 
coordonnées x; vers le système de coordonnée q' et inversement. 


En terme de distance infinitésimale dans l'espace, cette dernière que l'on appelle « métrique » est 
alors exprimée sous forme tensorielle de la manière suivante : 


ds? = lax} = "i a dq'dq’ = g ,dq'dq’ > g, = G Ok; 
dq' dq dq' dq 
ds’ =dq J" Jdq=  da=(dg',dg",dq )= g, =V 3 


En terme matricielle cela s'écrit également : 


Alors le tenseur métrique g; contravariant n'est ni plus ni moins que le produit de la transposée de 
la matrice Jacobienne par elle-même (uniquement dans le contexte d'un changement de 
coordonnées qui transforme l'espace tout entier en lui-même). On introduit dans ce contexte le 
tenseur métrique covariant gÏ comme l'inverse de la matrice correspondant au tenseur métrique 
contravariant : 


Lorsque les nouvelles coordonnées sont orthogonales alors le tenseur métrique est une matrice 
diagonale et la métrique s'exprime à l'aide des coefficients de lamé qui sont les termes diagonaux 


de ce tenseur: ds? =(h fla}  (nŸ =, = x a | 
dq' dq 


Dans le cas de coordonnées orthogonales alors le Laplacien s'écrit sous une forme assez simple : 


E E Er E FE z), 1 à 4 = 1 à (a r) 


Ve ôq'| h 6q'] Jg êg | h? ô) Jgop\ h og 


Dans le cas d'un système de coordonnées non-orthogonales, on peut utiliser l'expression plus 
générale du Laplacien de Beltrami, dans ce contexte de changement de variable : 


G=(e,)-3".3 > AT = a ô | TONE - 


| 
G) ôq q 
On note encore plus brièvement les dérivées partielles sous la forme : 


Ô 1 
= — = Ô. 
ôq k Videt(G) í | 


ô det(G) g'ar) 


C'est une forme extrêmement condensée mais en réalité si l'on devait l'exprimer avec les 
expressions générales de la transformation des coordonnées, cette expression est horriblement 
complexe, y compris dans le cadre d'un système de coordonnées orthogonales. 


Appliquons ces expressions au système de coordonnées non-orthogonale à 4 dimensions proposée 
par Pierre Humbert en 1936 dans son livre « Potentiels et prépotentiels » en page 34 : 


z= l-u’ -v F(p,0) 


Les fonctions F et ® s'entendent ici des fonctions similaires à celles assurant l'orthogonalité d'un 
système de coordonnées à 3 dimensions dont la paramétrisation est la suivante : 


oF(p,o)0F(p,0) | 0D(p,0) 0D(p,0) 
Ôp 06 Ôp 06 


=0 condition d'orthogonalité 


PC ET) 


S(p,0)= ea ET) 


ôo 


z=@(p,0) 


Calculons alors le Laplacien de Beltrami pour un tel changement de variable et pour cela 
commençons par les diverses matrices ou tenseurs métriques. C'est plus facile d'effectuer ces calcul 
avec Mathematica. La matrice Jacobienne s'écrit sous la forme : 


(p, ) 0 , (P,0) , F(,0) 
op 06 
J = P © 
,_Ft:9) Flo) pp ec) pur Eeo) 
1-u° -v l-u’ -v ô ôo 
0 0 0D(p,0) 0D(p,0) 
Ôp 06 
Le tenseur métrique contravariant est le suivant : 
ko) EE 5 o 
Sri) eo) l-u XElo.o)} 0 0 
l-u?’ -v l-u -v 
0 0 R(p,0) 0 
0 0 0  S(p,o) 


Le tenseur métrique covariant est le suivant (matrice inverse du tenseur métrique contravariant) : 


FE) Cey | 
uv : 1—v ; 0 0 
| Foo Foo) 
0 0 | 0 
R(P,6) | 
0 0 0 Ta 


Le déterminant du tenseur métrique contravariant est particulièrement simple : 


aele, )= ECD R(o,o)S(o,0)= faele, ] =- ELDE Jroa) 


On a donc tous les ingrédients pour calculer « facilement » la forme du Laplacien de Beltrami dans 
ce Cas : 


l-u’ OT 


1-u° y? u 


uv OT |. 
Vi-u?-v? ôv | 

uv oT 1-v* ôT | 
Ji- -v Ou re av J 


Sle,o) oT) © 2 
a) r) ôo (F(p.0) 


Rasha. | 
Ou 


V1-u? -v 
(F(p,0) R(P,6)S(p,0) 


el 


| 1 ô o) S(p,5)ôT |. 
| Fee Le (ee D EST) 


AT = 


R(p.5)5(p.0) 21 


| 1 Ô 
1-u°-v (p 
l-u? ƏT 


Vi-u?-v ôu 


V1-u° -v 


(Flp.o)} 


ô 
Ou 


uv OT 
i-u?-v? ôu 


2 (eoo Re) 


°c \ 


uv 


ôT 
V1-u° -v =) | 


L'équation de Laplace s'écrit donc : 


AT=0& 
es: 0 l-u? ƏT uv OT ô uv OT l-» ôT 
JR > S š 1 2 E -- H } 
l 7 (o 7 A u| J1-u?-v? ðu  J1-u?-v? ôv v| J1-u?-v? ðu  J1-u?-v? 0v 


=0 
Ô 


“fear enr) 2 


ôp 


(ero fen ar) 


Introduction des polynômes orthogonaux à bi-dimensionnels extension des polynômes de Legendre 
et de Chebyshev U 


Les polynômes de Legendre peuvent être définis par le développement en série de la fonction 
génératrice suivante : 


1 1 
J1-2a x+a° Ja xÝ +a(1-x°) n=0 


Les polynômes de Chebychev U de seconde espèce sont également définis par le développement en 
série de la fonction génératrice suivante : 


— = 
1—-2ax+a Pr 


Les polynômes de Legendre peuvent se construire à l'aide de la formule de dérivation suivante : 


P(x)=— 


| 2"n! 


ad” 
dx” 


(1) 


Les mathématiciens Charles Hermite et plus tard Paul Appell et J. Kampé de Fériet ont introduit 
puis travaillé sur des polynômes bidimensionnels et multidimensionnels orthogonaux, dans le cadre 
de l'étude des fonctions hypergéométriques à deux dimensions et plus. Une étude complète s'y 
trouve dans l'ouvrage de 1926 : « Fonctions hypergéométriques et hypersphériques, polynômes 
d'Hermite ». 


Ils ont introduit à ce propos les extensions suivantes à deux dimensions à partir de fonctions 
génératrice en x,y : 


1 m=+%0 n=+0 
= a"b'U,,,(x, y 
V(ax+b y 1P + (a? +b? À x? y?) 2 à #3) 


1 m=+0 n=+00 


D Da hs (x, y) 


1-2a x-2b y+ +b A Z 


Les polynômes Umn(X% y) sont l'extension des polynômes de Legendre, et les polynômes Vmn(X,y) 
sont l'extension des polynômes de Tchebychev de deuxième espèce. 


Les polynômes Umn(xX,y) peuvent se construire par une formule de dérivation : 


Ro 1 a | 2, : 2 Fe 
Un(x y)= 2 nm! dx" dy" (x | Yy 1) 


On peut aussi les construire selon les dérivations de la formule de Taylor appliquées à la fonction 
génératrice : 


1 qom 


1 
n!m! da” db” F ax—b y} + (a? p \i x? m 


U m.n (x, y) 


Les polynômes Vmn(X,y) peuvent se construire selon les dérivations de la formule de Taylor 
appliquées à la fonction génératrice : 


1 ge 1 
Fan (x, y) | | 


| n!m! da"db" |1-2a x-2b y+a° +b’ 


a=0 
b=0 


Les polynômes Umn(Xġ y) sont solutions du système suivant de deux équations aux dérivées 
partielles : 


ot (y) OU, (x, y) OÙ, (x, y) OU, (x, y) 
1 2 m,n ? m,n 2 2 m,n 2 1 m,n , 
ae Se a un | 
OU, , (x, y) OU, (x, y) OU, (x, y) OU, , (x, y) 
2 m,n > m,n ’ 2 m,n 2 1 m,n 0 
) ôy? i OxÔy (n y Ôy t i ôx (m 


nm H 1U „n (x, y)= 0 


n\n- I, »)= 0 


Les polynômes Vmn(X,y) sont solutions du système de deux équations aux dérivées partielles 
suivant : 


2 2 
(1 AE Fin (x, y) x ô Frs (x, y) ( 3) x OV nn (x, y) my OV nn (x, y) 


m (m+ n+ 2, (x, >)=0 


ox” á Ox0y ax ôy F 
2 2 
(1 >} CRE Ge y) xy ð ER (x, y) (m | 3) y OV n.n (x, y) n qV an (x, y) n (m Ln+ 2, , (x, y)= 0 
y Ox0y Ôy ox f 


Généralisation des polynômes de Gegenbauer à des polynômes bidimensionnels dû à Pierre 
Humbert 


Le mathématicien Pierre Humbert propose une généralisation tout à fait similaire des polynômes 
précédent, basée cette fois-ci sur la fonction génératrice des polynômes de Gegenbauer (attention 
les fonctions R ne sont pas polynômes) : 


s M=+0 n=+0 


(lax tb y—1) + (a° l p \i x? y) = X X a vue) 


m=0 n=0 


s-1 


(1 -x° — ye "© Farik, 5 


s—l 
Riy =l- -y 2 Us, (x) 


s 


(lax tb y-1) +4 (a° T x2 y) mO n0 
s+l m=+%0 n=+%0 


(1 2a x—2b y+a°+ CDI = >. Y'a"b' Vi, (x, y) 


m=0 n=0 


Cette généralisation est exposée dans l'ouvrage de 1926 des mathématiciens P.Appell et J.Kampé- 
de-Fériet « Fonctions hypergéométriques et hypersphériques, polynômes d'Hermite » dans un 
espace cartésien à un nombre quelconque de dimensions. 


Les polynômes {x y) peuvent se construire selon les dérivations de la formule de Taylor 


n+m s+l 
appliquées à la fonction génératrice : v} (x, y)= 1 A f 2a x-2b y+a° vY z | 
nm. aa 


Et de même pour les polynômes L,,,(x, y) et les fonctions R°,,,(x, y) associées : 


Un) ro [lex tb y—1) + (a i v (i x? »)i] 


“alm! da”db" 


Pons (CE E 2} 


n\m! da”db” 


On peut déduire les systèmes d'équations aux dérivées partielles de ces trois fonctions à partir 
d'une lecture attentive de l'ouvrage « Fonctions hypergéométriques et hypersphériques, polynômes 
d'Hermite » pour les fonctions R°,,(x,y) en page 277, formule (37), pour les polynômes U°mn(X,y) 
en page 278, formule (41) et pour les polynômes V°`mn(X%, y) en page 239, formules (32) et (32bis). 
Pour les polynômes (x, y), outre la réduction à deux dimensions, la subtilité consiste à remplacer 
dans la formule (32bis) page 239 le paramètre dimensionnel n par n+s-1, soit n par 1+s en deux 
dimensions, u étant la somme des indices du polynôme V,,,(x,y), ainsi qu'il est remarqué en page 
249, permettant d'obtenir les équations des V’,,.(x,y) en formule (53) page 250. 


Les systèmes d'équations aux dérivées partielles sont alors les suivants : 


ð y (x y) ð y; (x, y) V; (xy) V; (x, y) 
(x 2) a xX Go (n HS4 2)x à my d m (m+ n+l+ sV; (5 y)=0 

ð Vi (xy ô y; (x y) CARE ,y) av (x, y) S 
(1 y 2) xy að (m S+ 2)y d nx + n (m +n+1+ sV; (œ y)=0 

OR; ne }) © Ria (x (x 5) OR, lc y) OR; „(xX de y) OR, n (x, y) 
(i x a) 2 xy Go (n m+s 2)x T (m 1)| x à y > (m + 1) (m+ n+s IR: ne (x, y)= 0 

2 ps 2ps g g s 
(i y a)Ê R a% x,y) x Ô R; (x, y) m+s 2)y ôR (x.y) (n i|: ôR; (x,y) | y sei (n | 1) (m+ n+s 1); à ue y)= 0 
ôy OxÔy ôy ôx Ôy 
+ (x, y) 


(x fi A 


ox? 


(s fs (1- x? 


na (®©) 


ôx 


ôy 


nx ml x ur, is) =0 


U; y) (m 2) yra) (n tx Ur, (à, y). -( 


r OU: (xy) | G 2) Ura) (m 1 y e D). +(m tn\(m+ DU „œ » 


Ox0y 


OU: EU: 
| U; (x) a at, 


(x »{l FEU 


ôy 


2 


Ox0y 


(s f» (> 


jan») px OU, 109) (m y ni-y 2) 


ôy 


U* 


m,n 


0x 


A ne 


m + 


ne DU ls, »)}: 


Dans le cours du livre de Pierre Humbert « Potentiels et pré-potentiels », 1936, l'auteur utilise les 
fonctions R°,,{x,y) pour s=2 (introduit en page 36 et identifié en page ): 


R (xy)=(1-x- y RU, (a y) 


(x) 
(a x+b y—1) + (a? Ta x? y) 


m=+% n=+%0 


X J a”b" R, (xy) 


m=0 n=0 


Dont les équations aux dérivées partielles sont les suivantes : 


2 2 
m 1) £ OR, (x, y) | y OR, (x, n) 
ôx ôy 


2 2 
i|: OR,» (x, y) | y OR, (x, n) 
ôx Ôy 


z or 5) 


oR, n X, y ôR n X, y 
1 ED, en 


OxÔy ôx 


FR (x y) R? (x,y) oR? (x,y) 
1 Ran m,n m,n 
p-s E a PE) E gy 


RS 
= 


(+1) (m- n+ 1)R 2 (x y)= 0 


(n+1)(m+ n+ 1)R 2 (x,y)=0 


La sommation des deux équations aux dérivées partielles donne une nouvelle équations au 
dérivées partielles : 


PR y) 1 ER y) 

(Lx 7) ôx? ' ( y ) ôy? 
> (1 g r (1 >} R? (xy) 2x y 0 Re ») 2x 
X 


2p2 2 3 
Ta Rè (x,y) g y 2 Rama) 2 y © Run) + (m n+ 2)(m+ n+ 1)R PACS y)= 0 
x0y ôx ôy 


ô Riu (x, y) 2 ô Ri, (x, y) = ( 
y” Ox0y ôx ? ôy g 


m+n+1)(m+n+2) Rond) 


Séparabilité de l'équation de Laplace à 4 dimensions 


En revenant à l'équation de Laplace précédente : 


Ou Ne Ou Vie ôv) ðv \1-u?-v? Ou "M ôv 
ô 2 S(P,0 OT 0] 2 R(P,0 OT 
—| (F(p, H F\p, — |=0 
i eo Feer), 2 eoo er 


Développons les termes en u et v qui se simplifie et apparaît comme étant les termes de l'équation 
aux dérivées partielles de la fonction R°mn(X,y) : 


Jipoia] l-u? ôT w )-2{ uv ôT I-v 2) 


re ô l-u? ƏT uv oT) 0 uv oT 1-v? ôT 
Ou V1-u°—-y? ôu  J1-u?-v? Ov ðv V1-u?-v? ôu iwy ôv 
2 2 2 
= (1 7e) T+A }2T D RL DH 
Ou Ôv Ou Ov Ou Ov 


Si l'on sépare la solution recherchée de l'équation de Laplace, sous la forme : 
T(p,o,u,v)=R,, (u,v)Q(o,o) 


Alors la transformation de l'équation de Laplace donne une équation aux dérivées partielles pour la 
fonction Q s'écrit : 


G en iaaah, ô (Eloo) Heg calea E E T 


ap R(p,o)  ôp ôo S(p,o) ôo 


équation qui peut également s'écrire sous la forme : 


rleo) rloa) Bee] ale), moa) 2| rio.) [Eee] aleo) 
(m+n+1)(m+n+2)JR(b,0 )S(p,0 Q(p,0)+ 
©F(p,0) S(p,o) 0Q(p,0), 6F(p,0) p 


— F 
i op (aa) R(p,o) êp ôo Í Slo,o) ôo 


Comparativement elle est différente de l'équation aux dérivées partielles obtenus par la séparation 
en trois dimensions dans un système axisymétrique : 


= F0) F0.) PE lea), Foo) Erloa) PET alpa ANESSE 


Extrait du calcul de René Lagrange pour le système axisymétrique 3D 


Nous allons maintenant analyser l'équation du scalaire U d'une manière similaire à l'analyse 
menée par un article assez fondamentale de Renée Lagrange de 1939 :« Les familles de surfaces 
de révolution qui possèdent des harmoniques », dans lequel l'auteur analyse les conditions pour 
lesquelles la séparabilité simple ou la R-séparabilité dans un système axisymétrique est rendu 
possible et l'équation différentielle ordinaire qui sous tend une fonction de révolution à la base de 
la définition du système de coordonnées orthogonales dans le plan méridien. 


De la condition d'orthogonalité : OF OF _ OP êP =0 on peut déduire l'existence d'une fonction 


ôp ôo ôp ôc 
AÀ(p,o) telle que : 
2 2 
op or ee r(p,o) =E) +x(2) 
op ôp 2 ôp ôo ap 00 R 
£ =- -= Alp,o)> > 5 PEEN 
Cr ôD __10r 1EY EF) ŸS 
ôo ôo oo aap [Soala al) 


2 2 
2 oF oF 1 
o Hy)’ = Zi VRS = 
n note ( 1) (2) (2) > A4) 


Et l'équation au dérivées partielles dans le plan méridien s'écrit alors : 


(re) {2 (n+1ÿ 2-0 


ôp\ À ôp) ôs ôs 


La séparabilité est assurée si la fonction prend la forme Q(p,0)= f(p,0)N(p)T (o), où f est une 
fonction de séparation, et si l'équation se sépare en deux équations différentielle ordinaire. 
D'abord l'équation précédente s'écrit sous la forme : 


F [z ÔN ô [ue LE) | EE ôT ô (osla Fr} Ll2 (324) ô 6) (m+ i 
AN | õp ôp êp À T | ðo% ôo ðo f\ðp\ à õp) ĉc o FAH, 


ag 1 [zr | ON ð [ue LE) | SE | ôT ô oela F) 1 k (2%); ð Le) Get : 
AN | ð?  Ôp ôp À T | ôo? ðc ðc fF|õp\ à p) ĉc ô F’AH, 


La première condition de séparabilité est que les deux expressions soit la somme de termes, l'un 
dépendant uniquement de p et l'autre de o : 


zou LE )-dto) blo) aar eoad cerner) 


0 


En prenant une forme de la fonction À telle que : to ), alors le système 
d'équation pour les fonctions s'écrit : 


0D_n(p)oF 1 ôb 1 ©F 
p  t(o)ôc n(p) êp t(o)ôo 
D (o) ôF 1 ôb _ 1 ôF 
ôo n(p)@ (o)6o  nlp)ôp 


En prenant des primitives des fonctions n et t comme variable à la place de p et © alors le système 
s'écrit : 


Système qui équivaut à choisir une fonction À=1 par ce changement de variable. Donc oin n 
restreint pas la généralité de l'étude en posant : À=1. 


Si l'on pose À=1, la conséquence immédiate est que les deux fonctions F et ® sont toutes deux 
harmoniques (et conjuguées l'une de l'autre par la condition d'orthogonalité) : 
ô ôF PP PD 
õp ôo a ôp? ôo’ 
ôb © @F F 
ôs êp (ðP ôo’ 


0 


Remarque : c'est l'argument classique qui permet de choisir comme système de coordonnées 
orthogonales, parmi les parties réelles et imaginaires des fonctions analytiques du plan complexe. 


Les conditions pour la séparabilité sont Log(Ff*)= n(e)+a(o). par ailleurs on peut noter avec 
l'hypothèse À=1 (ne restreint pas la généralité du propos selon R.Lagrange) que : 


R(p.6)=5(p,0)= A, =(H,)° (2) | (£) 


Ôp 06 


Ainsi la deuxième condition nécessaire pour que l'équation soit séparable est que le terme : 


1 fr) 3 de) (m+1)?A, 
Ff |ôp\ p) 60o\ ôs F? 


soit également somme d'un terme dépendant uniquement de p et d'un terme uniquement de o, ce 
qui conduit si cela doit être vrai quelque soit m à la double condition : 


FE G 21) | Z G ZL)| = slp) ste) 


Â = g(p)+h (o) 


La première de ces deux conditions est immédiatement vérifiée, en posant d'abord 
g l r o s(o 
p.(p)=2[dp r(p) a(s)=2[dc slo)=> f F (p)Qlde slo) 


Puis en écrivant: 


of fF, _ fôlog K a | 
o Meg HN € SOFIANE 


of fF _ J0Log “0: J 
ðo 2F ôo HEIE 2 ôo ro) * ôo (WF AN 
Cf. fõlogF f 2e) Mr EE roA S D'Logr _ f f an 1 H 1 °F 


ôp 2 ô sf ôp ôp op  @plFèôp) Pla) Fp 


ôp\ ôp ôp” 2| êp ôp 2 ôp” 4l ôp 4F°\ õp) 2F8p 
o( of 1 (FF 1 0F A 

F =F +s'(o )+ 
2 2) (F7 2F do so) Ce) 


Ce qui donne : 


A) E E rene secte 


Cela démontre donc la première décomposition. 


Si nous revenons à notre équation équation aux dérivées partielles dans le contexte de la 
séparabilité à 4 dimensions : 


Alors l'équation de la fonction Q s'écrit : 


z ner tie) Le HF.) ape) 222) (m+n+1)(m+n+2) A(p.0) = 0 


p| A(p,0) ôp do 06 A(e,6X4,(0,0)) 


Que l'on va simplifier en omettant les arguments dans les fonctions, soit : 


2 
CA REA (ra) (m+n+1)(m+n+2) i 
dp\ À Ôp) ôo Glo 


(rot) sE 2 1 ô’ LogF [E] voeon >a l er) l 00) 


Si la fonction Q est définit au moyen d'une formule de séparation : Q(p,0)= f(p,0)N(p]T(o), où f 
est une fonction de séparation, alors l'équation précédente s'écrit sous la forme : 


F? [zx ÔN à pee) | mafon ƏN ô tolr) 


AN | ôp? ôp ôp 4 T | ôo ôo ôs 
Lo Pf) 2 (r2) (m+n+1)(m+n+2)_ o 
 flõp\ à p) 60 ðc 2H? 
2 2 r2 2, 
1 CN ON te E „2 CEREN a 
AN | 0p°?  ôp ôp À T| ðo? ôo ðc 
1 jaf Poff) (raf) (n+n+1)(m+n+2)_ o 
 F°flõp| à ôp) 60 ðc AHP F’? 


Pour que cette dernière équation se décompose en deux équations différentielles ordinaires selon p 
et o, il est nécessaire que : 


Zro 
CE ] indépendante de o 


op À 
L (Logla F’ f? ) indépendante de p 
o 


en adoptant le même raisonnement que Renée Lagrange, page 286, et son point de vue sans 
restriction de généralité en posant À=1, il vient que l'expression F°F doit être un produit d'une 


fonction de p et d'une fonction de o, soit donc : Log(F’ f° )= pi(e)+a; (e), Par ailleurs on peut noter 
avec l'hypothèse À=1 que : 


R(p,0)=S(p,0)= A, =(H,)° (£) (#) 


op) ôo 


Ainsi la deuxième condition nécessaire pour que l'équation soit séparable est que le terme : 


1 fra) 3 (32) (m+n+1)(m+n+2)A, 


F?flôpl p] ôs ôs F? 


soit également somme d'un terme dépendant uniquement de p et d'un terme uniquement de o, ce 
qui conduit si cela doit être vrai quelque soit m et n à la double condition : 


ne eoe sono) 


2, = g(p)+h(o) 


Tout comme dans l'article la première de ces est la conséquence de deux 


P(p)=2{dp r(p) c)=2[do s(o j= f=e fdp ro) oldo slo) 


décompositions. En posant : . Ici si sous 


nous posons F = VG, nous avons : 


CC) e ET eeo 

Soit : 

F7 ale G (ré) PE COS EN) oeo PTT) 
er) (28) Ja) +28) je 


(EEE) GE] rotto este 


õp? ôo?) 4F°|\6p °0 


PER 
~ 
gpa, 
D 
Ner 
Ne 
N 
a, 
PER 
Q 
Ner 
e; 
Ca 
PES 
Q 
Nr 
Ner 
N 


Ce qui démontre donc la première décomposition. 


Pour qu'il y ait séparabilité, il est donc nécessaire que les deux conditions suivantes soit remplies : 


Log(r° f°)= r(e)+a.(o) 
AO (p)+h(0) 


Les intégrales primitives de la fonction F, qui par ailleurs est une fonction harmonique de p et o 
sont donc telles que : 


AF=0 
CON OC Or 


Il y a donc deux résultats : 

- la « bonne nouvelle » c'est que l'on peut construire un système de coordonnées non-orthogonales 

à 4 dimensions selon la paramétrisation précédente pour lequel l'équation de Laplace est 

séparable. Pour cette dernière on utilise l'extension de polynômes de Legendre à deux dimensions 

— la « mauvaise nouvelle » c'est que les systèmes de coordonnées sur l'équivalent d'un plan 

méridien à 4 dimensions sont exactement les mêmes que ceux sur d'un système 
axisymétrique sur son plan méridien : cartesien, polaire, elliptique, parabolique, et leur 
inversion par rapport à l'origine des coordonnées, courbes quartiques de la famille des 
cyclides de Dupin (ce sont des sections de cyclides de révolution, René Lagrange les 
désignent par « courbes cartésiennes », Limaçon de Pascal, cardioide, courbes à trois foyers 
, …). Dommage que l'étude ne puissent introduire de nouvelles catégories de courbes 


à finir : introduire la forme de la fonction de révolution et de séparation, voir travaux de René 
Lagrange 


